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Abstract 

The  flow  distribution  features  in  U-type  layers  or  stacks  with  certain  practical  configurations  have  been  investigated  analytically.  The  formulations 
suggest  general  designing  strategies  to  improve  the  flow  uniformity,  and  narrowing  some  positions  of  the  channels  proves  effective.  The  flow 
uniformity  among  the  layers  in  a  U-type  stack  is  relatively  easy  to  achieve  in  comparison  with  that  among  the  channels  in  a  U-type  layer,  due  to 
the  large  stack  headers  and  low-pressure  loss  in  them.  CFD  simulations  confirm  the  formulations,  and  the  discrepancies  between  the  analytical  and 
CFD  results  have  been  attributed  to  the  ignored  factors  during  the  analytical  formulations. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  flow-distribution  uniformity  is  a  basic  target  in  fuel¬ 
cell  design  and  manipulations,  and  various  configurations 
have  emerged  to  improve  the  performance.  The  most  popular 
method  for  design  optimization  is  usually  based  on  numeri¬ 
cal  techniques,  adopting  commercial  CFD  (computational  fluid 
dynamics)  codes  [1-4],  for  the  cases  are  normally  too  com¬ 
plicated  to  treat  analytically.  However,  analytical  solutions  are 
extremely  desirable  if  available  because  of  their  straightfor¬ 
ward  and  fast  solving  and  possibly  revealing  definite  relations 
among  quantities.  Fortunately,  certain  configurations  have  been 
reported  with  analytical  solutions,  based  on  some  simplifica¬ 
tions,  including  the  well-known  Z-type  and  U-type  parallel 
channels  [5,6].  Based  on  analytical  formulations,  Sung  [7]  has 
further  proposed  optimized  header  shapes  to  improve  the  flow 
distributions. 

Thanks  to  the  difficulties  in  treating  complex  pressure  losses, 
the  analytical  formulations  reported  so  far  are  mainly  empha¬ 
sized  on  relatively  simple  configurations,  and  at  the  layer  level. 
Meanwhile,  more  complex  configurations  are  popular  in  prac¬ 
tice,  and  analysis  on  a  fuel-cell  stack  is  commonly  as  necessary 
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as  on  a  layer,  for  flow  distribution  among  the  layers  in  a  large 
stack  can  play  a  controlling  role  on  the  overall  performance. 
With  certain  simplifications  and  assumptions,  the  present  work 
attempts  to  formulate  analytically  in  some  more  complex  cases. 
Since  the  U-type  configuration  is  more  efficient  than  the  Z-type 
one  in  improving  the  flow  uniformity  [5],  and  more  widely 
adopted  in  fuel-cell  stacks,  U-type  layers  and  stacks  are  the 
focus  here.  Comparisons  with  CFD  simulations  will  also  be 
discussed. 

2.  Formulation 

Certain  assumptions  are  employed  in  the  following  formu¬ 
lations.  The  density  p  and  dynamic  viscosity  [i  of  the  fluid  are 
assumed  to  be  constant,  i.e.,  the  flow  is  assumed  to  be  incom¬ 
pressible  and  the  temperature  is  assumed  to  be  uniform.  The 
gas  consumption,  leakage  or  transformation  is  not  taken  into 
account,  and  the  mass  conservation  is  assumed  in  the  whole  flow 
field.  The  flow  is  assumed  to  be  laminar,  and  the  inertial  term 
is  neglected  during  considering  the  momentum  balance.  The 
local  losses  are  neglected  in  all  cases,  including  the  dividing, 
combining,  expanding  and  contracting  ones.  Especially  in  the 
stack  cases,  the  pressure  loss  in  the  layer  headers  has  also  been 
neglected  to  simplify  the  formulation.  And  for  convenience,  the 
feed  header  in  present  work  is  assumed  to  have  the  same  size  as 
the  exhaust  header. 
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Nomenclature 

a  aspect  ratio 

A  cross-sectional  area 

D  hydraulic  diameter 

/  friction  factor 

F  flow  parameter  index 

g(F )  function  of  F  in  Eq.  (40),  etc. 

G\  constant  in  Eq.  (14a),  etc. 

G2  constant  in  Eq.  (14b),  etc. 

H  height  or  depth 

K()  intermediate  quantity  in  Eq.  (6),  etc. 

K\  dimensionless  quantity  in  Eq.  (9a),  etc. 

K2  dimensionless  quantity  in  Eq.  (9b),  etc. 

/  overall  loss  coefficient 

L  length 

N  number 

p  pressure 

P  perimeter 

r  constant  in  Eq.  (15),  etc. 

Re  Reynolds  number 

s  constant  in  Eq.  (15),  etc. 

u  velocity 

W  width 

x  Cartesian  coordinate 

Greek  letters 

Afric  frictional  loss  coefficient 

/x  dynamic  viscosity 

p  density 

ft  turning  loss  factor 

Subscripts 
avg  average 

c  channel 

cl  channel  parti 

c2  channel  part  2 

e  exhaust 

f  feed 

h  header 

in  inlet 

layer  layer 

max  maximum 

min  minimum 

r  rib 

stack  stack 

Cap 

dimensionless 


Mass  balance  requires 


^h-layer 


dl/f_iayer 


Nc 


dVf_iayer  ^h-layer 

for  the  feed  header  and 


UrAc 


^h-layer 


dwe-i 


e-layer 


Nc 


dx 


e-layer 


^h-layer 


UqAq 


(la) 


(lb) 


for  the  exhaust  header,  where  A,  u ,  x,  N  and  L  stand  for  the 
cross-sectional  area,  the  velocity,  the  distance  from  the  inlet,  the 
number,  and  the  length,  respectively,  and  the  subscripts  “h”,  “f”, 
“e”,  “layer”  and  “c”  means  “header”,  “feed”,  “exhaust”,  “layer” 
and  “channel”,  respectively.  With  the  ignorance  of  the  inertial 
term  [5],  the  momentum  balance  requires 


“(P^f-layer)  —  “ 


d  Pf-layer  Ai-layer  /h-layer  Puf -layer 


dXf_iayer  fi-^f-layer 

for  the  feed  header  and 


2  A  h-layer 


(lc) 


dx 


e-layer 


(Pu  e-layer)  ~  _ 


d  Pe-layer  ^h-layer  /h-layer  P^g-layer 


dx{ 


+ 


e-layer 


2Ah_iayer 


(id) 


2.7.  U-type  layer  with  uniform  channels 

The  formulation  in  this  section  involves  the  U-type  config¬ 
uration  containing  channels  with  uniform  shape  (Fig.  1,  called 
“casel-0”),  basically  according  to  the  work  in  Ref.  [5]. 


for  the  exhaust  header,  where  p,  P  and /represent  the  pressure, 
the  perimeter  and  the  friction  factor,  respectively.  Considering 

1  2 

Pf-layer  —  Pe-layer  —  Apf-e-layer  —  ^c^P^c  (2) 
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while 

k  =  £t-f  H”  £t-e  H”  ^fric-c  ~  ^fric-c 

4/cLc 


7-fric-c  — 


Dc 


(3) 

(4) 


where  lc  is  the  overall  loss  coefficient  in  the  channel,  £t-f  and 
f  t-e  are  the  turning  loss  factors  for  flow  dividing  and  combining, 
Afric-c  is  the  frictional  loss  coefficient  in  the  channel,  and  D  is  the 
hydraulic  diameter.  Neglecting  ft_f  and  ft_e  in  Eq.  (3)  is  based 
on  the  low  Reynolds  number  (laminar  flow)  and  the  high  length 
to  diameter  ratio  of  the  channel  (. Lc/Dc  >  100)  [5]  investigated 
here.  Using  the  empirical  correlation 


Ref  =  13.84  +  10.38  exp 


-3.4 


(5) 


where  a  is  the  channel  aspect  ratio.  According  to  Eqs.  (2)-(4) 
and  the  definition  of  Reynolds  number,  the  relation  between  uc 
and  A/?f_e_iayer  can  be  derived  as 


A/?f-e-layer  —  1 


(2tiLc(Ref)c 

V  Do 2 


uc 

K~0 


(6) 


where  fi  is  the  dynamic  viscosity  of  the  fluid.  With  the  following 
nondimensionalizaion: 

(7a) 
(7b) 

P  =  -  (7c) 


T^h-layer 

U 

^in-layer 

P 


PU  in-layer 

Eqs.  (la)-(ld)  evolve  as  a  set  of  ordinary  differential  equations: 

d/if layer 


dx 

dw  e-layer 

dx 


—  K\  A ^f-e-layer 
=  ~K\A  ^f-e-layer 


du  f-lay er  d/^f-layer  v  ^  A 

2l/f_iayer  T7  I  ~T7  I”  A2Wf-layer  —  0 


2  fi 


e-layer 


where 


dx  dx 

dw  e-layer  d/?e_iayer 


dx 


+ 


dx 


K2  fi  e-layer  —  0 


m 

(8b) 

(8c) 

(8d) 


Ncpu  in-layer Ac  NcpU  in-layer 

K\  =  - - - Kq  = 


^h-layer 

ACD2C 


2  fi 


(Re  f) cT'c^h-layer 


m 


k2  = 


T^h-layer  p(R  £/)h-layer^h-layer 
2Ah-layer  ^h-layer  PP  in-layer 
2/X  (^^/^)h-layer^h-layer 


Pu  in-layer 


D\ 


(9b) 


According  to  the  continuity  equation 

fif-layer  =  fi  e-layer  (Id) 

A  second-order  differential  equation  is  derived  as  following 


d  U  f_iay er 

dx2 


2^1^2^f-layer  —  0 


Using  the  boundary  conditions 

X  =  0,  A  f-lay  er  =  1 
X  =  1 ,  fif- layer  =  0 

the  resultant  solution  is 
fif-layer  =  G\e  "  G2C 
where 


G2 


and 


es  —  er 
er 

er  —  es 


(11) 

(12a.) 

(12b) 

(13) 

(14a) 

(14b) 


r  =  -s=y/2KiK2  (15) 

Meanwhile,  the  following  results  can  be  obtained 
G\rerx  +  G2ses* 


^Pf-e-layer  — 


mc  =  -- 


Ki 

G\rerx  +  G2sesx 

~N* 


(16) 


(17) 


where  m  is  the  dimensionless  relative  flow  rate.  The  follow¬ 
ing  parameter  can  be  utilized  to  evaluate  the  uniformity  of  flow 
distribution 


Fc  = 


Pic-  max  Pic-  r 


TYlr-  r 


(18) 


2.2.  U-type  layer  with  complex  channels 


h-layer 


During  practical  manufacturing,  complex  channels  might  be 
selected  owing  to  the  technical  limitations  or  cost  considerations. 
If  certain  shape  change  can  improve  the  flow  distribution,  it  is 
certainly  further  desirable. 

The  first  case  considered  here  is  that  each  channel  consists  of 
three  parts  (Fig.  2,  called  “easel- 1”).  The  two  parts  connected 
respectively  to  the  feed  and  exhaust  headers  are  of  the  same  size 
and  called  “cl”  while  the  middle  channel  is  called  “c2”.  Hence, 
this  kind  of  channel  is  named  as  “c  1  +  c2  +  c  1”  type  in  the  present 
work.  For  convenience,  the  “cl”  here  and  hereafter  has  the  same 
width  as  “c2”. 

When  the  local  losses  described  in  Eq.  (3)  and  those 
across  the  two  channel  parts  (expanding  and  contracting 
losses)  are  neglected,  the  following  relationship  can  be 
derived: 
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Fig.  2.  The  U-type  layer  with  “cl  +  c2  +  cl”  channels  (casel-1). 


Apf-e-layer  —  2  [A 


( Lcl(Ref)cl  ,  Lc2(Ref)c2  Aci 

V  Dcl 


+ 


D 


AC2 


Uq\ 


Uc\ 

K0 


(19) 


where  Lc\  =LC—  Lc 2,  is  the  total  length  of  “cl”  in  a  channel.  In 
this  case,  mass  balance  gives  the  following  relations: 


Ah-layer 

Ah-layer 


d^f-iayer 

cL^f-layer 

dw  e-layer 


Nc 


f^h-layer 

Nc 


wcl  AC1 
*Al  AC1 


(20a) 


(20b) 


cL^e-layer  f^h-layer 

while  the  momentum  balance  provides  the  same  equations  as 
Eqs.  (lc)  and  (Id).  Eqs.  (8a)-(8d)  can  also  be  derived  with 

^  in-lay  erAcl  ^  Ncpu in-layer 

K\  =  - - - Aq  =  - 


Ah-layer 


2/X 


KLci(Re  /)ci/ Ac\Dcl  +  Lc2(Re  f)c2/ Ac2Dc2J  Ah-layer 

(21) 


and  the  same  expression  of  K2  in  Eq.  (9b).  After  solving,  the 
same  results  as  expressed  in  Eqs.  (16)  and  (17)  can  be  obtained. 

Compared  with  case  1-0,  it  is  observed  that  partitioning  the 
channel  into  different  parts  has  brought  extra  term  into  the 
pressure  loss  along  the  channels,  helping  to  tailor  the  flow  dis¬ 
tribution. 

The  next  case  for  investigation  is  that  the  channel  consists  of 
only  two  parts  (“cl”  and  “c2”)  connecting  the  inlet  and  exhaust 


header,  respectively  (Fig.  3,  called  “case  1-2”),  and  such  channel 
is  called  “cl  +  c2”  type  hereinafter. 

In  this  case,  the  same  relationship  as  Eq.  (19)  holds  if  the 
local  losses  are  neglected.  The  mass  balance  gives  Eq.  (20a)  for 
the  feed  header,  and  the  following 

di/e.iayer  Nc  Nc 

Ah-layer“j  —  ~~  uclAc2  —  ~~  UC\AC\  (22) 

UTe-iayer  ^h-layer  ^h-layer 

for  the  exhaust  header  due  to  uc\  Ac\  =  uc 2  AC2.The  further  for¬ 
mulation  is  the  same  as  the  above  casel-1  with  K\  expressed  in 
Eq.  (21),  K2  in  Eq.  (9b),  and  the  results  in  Eqs.  (16)  and  (17). 

More  complex  channel  shape  can  also  be  analyzed  using  the 
same  formulating  procedure  as  above.  For  any  types  of  chan¬ 
nels  containing  different  parts,  the  frictional  loss  corresponding 
to  each  part  can  be  considered,  respectively,  and  the  only  cor¬ 
rection  is  within  the  K\  expression  when  the  turning,  dividing, 
combining,  expanding  and  contracting  losses  are  neglected. 

2.3.  U-type  stack  with  parallel  layers 

U-type  stacks  are  frequently  encountered,  for  the  feeding  and 
exhausting  are  usually  conducted  on  the  same  manifold  plat¬ 
form  to  easy  the  gas  manipulation.  The  flow  distribution  among 
layers  is  never  less  important  than  that  among  channels.  Fortu¬ 
nately,  U-type  stacks  are  very  analogous  to  U-type  layers,  and  the 
formulation  becomes  timesaving  on  the  basis  of  the  preceding 
sections. 

First,  uniform  channels  are  investigated.  A  fuel-cell  stack  is 
considered  with  layers  forming  a  U-type  configuration,  while 
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each  layer  contains  the  same  uniform  channels  (Fig.  4,  called 
“case2-0”). 

In  order  to  investigate  the  flow  distribution  among  the  lay¬ 
ers,  the  mass  balance  and  momentum  balance  in  the  stack  are 
considered  as  follows 

A  di/f_stack  _  flayer  ^  A 

^h-stack  “j  —  flayer  Mayer 

QXf-  stack  L  h-  stack 


y—^—uc-avgNcAc  (23a) 

f-'h-  stack 


and  A/?f_e_stack  can  be  derived  as 


1  2 

Apf-e-stack  —  hwy&r~Pu  layer 
_  M^c-avg 

“  Ko 

based  on  the  assumption  that 
4/cLc 


MMc-avg 


2/iLc(Ref)c 

Dl 


(24) 


Aayer  ~  Ifric-layer  —  Me 


Dc 


(25) 


and  wiayer  =  ^c-avg  •  After  nondimensionalization  as  follows, 
x 


x  = 


u  = 


Ph- stack 
U 

U  in-  stack 

P 

P  =  - 

PU  in- stack 

the  following  equations  can  be  achieved. 

dllf-  stack 


dl 

duQ- stack 
dl 


P/ A  pf_e_stack 
=  —  K\  Apf-e-stack 


d/if-stack  dpf_stack  .  „  ^  ,, 

ZU f_stack - — - 1 - — - V  A2Wf_stack  =  0 


2  We-' 


e- stack 


where 


dl  dl 

du  e-stack  .  dpe_s  tack 


dl 


+ 


dl 


^2  We- stack  —  0 


^  MayerP^in-stackAc  Mayer  Pw  in- stack 

Al  =  - - - Kq  =  - 


Ah- stack 

AcD2c 


2/x 


(Re  f  )c  Lc  A  h-stack 


(26a) 

(26b) 

(26c) 

(27a) 

(27b) 

(27c) 

(27d) 


(28a) 


A  dw  e-stack 

^h-stack 

^e-stack 


Mayer 
P  h-stack 
Mayer 
P  h-stack 


^  layer  A  layer 
wc-avgM  Ac 


(23b) 


Avf_stack 


(2  \  _  ^Pf-stack 

IP^f-stack)  —  ~~aZ 

dVf_  stack 


Ph- stack  /h-stack PU f- stack 

2  A  h-stack 


(23c) 


dve-: 


e- stack 


(p^e-stac 


k)  = 


4 Pe- stack  Ph- stack  /h-stack  Pu  e_ stack 


dve-< 


e-stack 


2Ah-stack 


(23d) 


where  the  subscript  “stack”  means  the  stack,  Mayer  is  the  number 
of  layers,  Nc  is  the  number  of  channels  in  one  layer,  wiayer  is 
considered  as  wc-avg  (the  average  uc ,  wc-avg  =  (l/M)S^iwc-/)? 
and  A  layer =  Mc  Ac.  Meanwhile,  the  relationship  between  wc-avg 


„  Ph- stack  P' (  P^Z/h- stack  P  h-stack 

a  2  =  - 

2  A  h-  stack  P^h-  stack  PU  in-  stack 

(  P  g/)h- stack  P  h-stack 

n 2 

^h-stack 

Using  the  continuity  condition 


uf- stack  —  we- stack 


and  the  boundary  conditions 

1  =  0,  Wf- stack  =  1 
1=1,  stack  =  0 

the  solution  is 

wf- stack  =  G  \e  +  G^e 

and 

G\rer*  +  G2ses* 

A  Pf-e-stack  —  ~ 


2/1 

Pu  in- stack 


(28b) 

(29) 

(30a) 

(30b) 

(31) 

(32) 
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Fig.  5.  The  U-type  stack  with  “cl  +  c2  +  cl”  channels  (case2-l). 


^  layer  — 


G\rerx  +  G2sesx 
Mayer 


(33) 


with  the  same  definitions  as  in  Eqs.  (14)  and  (15),  and  the  flow 
distribution  among  the  layers  can  be  evaluated  using 


^layer  — 


PI  layer-  max  M  layer-  min 

PI  layer-  max 


(34) 


Given  a  glance  at  the  expressions  of  K\  and  K 2  in  Eqs.  (28a) 
and  (28b),  it  is  easy  to  find  that  besides  channels,  the  flow  dis¬ 
tribution  among  layers  is  dominated  by  the  stack  headers  here 
instead  of  layer  headers  in  the  layer  cases.  This  makes  it  easier 
to  achieve  a  better  flow  distribution  among  layers  than  among 
channels  due  to  the  usually  larger  size  of  stack  headers  than  layer 
headers. 

Next,  we  begin  to  analyze  the  stacks  containing  complex 
channels.  The  first  case  is  that  each  channel  consists  of  three 
parts,  i.e.,  “cl  +c2  +  cl”  (Fig.  5,  called  “case2-l”). 

Similar  to  the  above  analyses,  the  following  equations  can  be 
deduced  based  on  mass  balance: 


Ah- stack 


fl^f- stack 
dTf_stack 


Mayer 
^h- stack 
Mayer 
^h- stack 


w  layer  ^  layer 
wcl-avgMAci 


(35a) 


A  dwe_stack  _  Mayer  ^  A 

^h- stack  ~  —  ~  ^  layer  layer 

d*e- stack  Eh- stack 

=  --^L«cl.avgJVcAcl  (35b) 

^h- stack 

where  zqayer  is  considered  as  wcl-avg  (the  average  uc  1,  wcl-avg  = 
(l / NC)Y^^=\U d -i),  and  Aiayer  =  MAci .  The  momentum  balance 
gives  the  same  equations  as  Eqs.  (23c)  and  (23d).  Once  only  the 


frictional  loss  is  considered,  the  relationship  between  wcl-avg  and 
A/7f-e- stack  can  be  derived  as 


M?f-e-stack  —  2/X 


(Lcl(Ref)cl  Lc2(Ref)c2  Aci  \ 
V  Dci2  Dc22  Ac2) 


M^cl-avg 


M^cl-avg 

K0 


(36) 


Subsequently,  the  deduced  differential  equations  are  same  as 
Eqs.  (27a)-(27d)  with  nondimensionalization  as  expressed  in 
Eqs.  (26a)-(26c)  and  K2  in  Eq.  (28b),  while  K\  as  follows. 


v  Mayer  Pu  in-stack  Aci  Mayer  PP  in- stack 

K\  =  - - - Aq  =  - 


Ah- stack 


2/X 


\Lc\(Re  /)ci/ Ac\Dcl  +  Lc2(Re  f)c2/ Ac2Dc2J  Ah-Stack 

(37) 


Based  on  the  continuity  and  boundary  conditions  shown  in  Eqs. 
(30),  the  solutions  give  the  same  expressions  as  in  Eqs.  (32) 
and  (33). 

Secondly,  each  channel  is  considered  to  be  the  “cl  +  c2”  type 
(Fig.  6,  called  “case2-2”),  with  “c2”  connecting  the  exhaust 
header. 

The  mass  balance  for  the  exhaust  header  gives 

dwe_ stack  Mayer  AT  . 

^h-stack“j  —  —  ~  wc2-avgM^c2 

CVe-stack  ^h-stack 

=  -2^L«ci-avg^cAci  (38) 

^h- stack 

according  to  wci_avg  Aci  =  wC2-avgAC2,  where  uc2  = 
(1/M)Z)£imc2-j-  Further  formulation  is  the  same  as  the 
above  case2-l. 
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2.4.  U-type  stack  with  U-type  layers 

Now,  we  explore  a  more  complex  configuration:  a  U-type 
stack  built  by  U-type  layers.  Firstly,  the  channels  are  considered 
to  be  uniform  (Fig.  7,  called  “case3-0”).  The  same  equations  as 
in  Eqs.  (23a)-(23d)  can  be  obtained  when  assuming  the  wiayer  as 

Nc 

wc_av g  =  (l/Nc)Yu»j.  And  when  neglecting  the  pressure  loss 

i—\ 

in  the  layer  headers,  the  same  expression  of  K\  as  Eq.  (28a)  can 
also  be  deduced.  With  the  same  expression  of  K2  as  Eq.  (28b), 
the  results  are  the  same  as  Eqs.  (32)  and  (33). 

When  the  channels  are  of  the  “cl  +  c2  +  cl”  type  (Fig.  8, 
called  “case3-l”),  the  formulation  is  same  as  “case2-l”.  And  if 
the  channels  belong  to  the  “cl  +  c2”  type  (Fig.  9,  called  “case3- 
2”),  the  formulation  is  same  as  in  “case2-2”. 


3.  Discussion  and  validation 


3.1.  Dependence  of  F  on  K\  and  K2 


As  a  simplification,  mmax  and  mmin  may  be  assigned  to  the 
values  of  m  at  x  =  0  and  1,  respectively,  during  the  calculation 
of  F.  Based  on  the  above  simplification  and  Eqs.  (14),  (15),  (17) 
and  (33), 


^  t  G\rer  +  G2ses  1  2er 

F  =  1 - =  1 - Tf- 

Gir  +  G2s  1  +  e2r 

Hence, 


r 


1  +  vn  2  -  f) 


=  g(F) 


(39) 


And  according  to  Eq.  ( 1 5),  the  following  relation  can  be  deduced, 
log  K2  +  log  Ki  =  2  log  (In  g(F))  -  log  2  (41) 

For  a  given  F ,  log  K2  ~  log  K\  appears  as  a  line  with  slope  of 

—  1,  as  shown  in  Fig.  10.  In  other  words,  a  line  with  slope  of 

—  1  in  the  log  K2  ~  log  K\  plot  represents  a  same  F  value.  More 
accurately,  mmax  and  mmax  can  be  calculated  as  the  mean  values 
in  the  first  and  last  channels  (or  layers),  respectively,  for  layer 
cases  (or  stack  cases).  However,  for  a  definite  configuration,  the 


Lc\ 


Lc2 


Fig.  9.  The  U-type  stack  with  U-type  layers  with  “cl  +c2”  channels  (case3-2). 


e' 


1  -  F 


(40) 
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Fig.  10.  Dependence  of  F  on  K\  and  K2,  based  on  simplified  calculations. 


3.3.  Applications  and  validations 

Some  examples  of  the  above-mentioned  cases  are  given  in 
Table  1,  along  with  the  corresponding  analytical  results.  For 
convenience,  the  F  values  in  Table  1  were  calculated  using  the 
mass  flow  rates  at  the  channel  centers.  The  results  are  in  accor¬ 
dance  with  the  discussions  in  the  above  section,  and  it  is  obvious 
that  changes  in  the  channel  size  can  tailor  the  flow  uniformity 
effectively. 

The  validation  of  the  analytical  formulations  depends  on  the 
correctness  of  the  related  assumptions.  During  the  formulation, 
several  assumptions  have  been  made,  analyses  of  which  will  aid 
to  understand  the  corresponding  validations. 

An  assumption  is  the  neglecting  of  the  inertial  term  during 
building  the  momentum  balance  equations.  In  the  case  of  tur¬ 
bulent  flow,  Bassiouny  and  Martin  [8,9]  neglected  the  frictional 


relation  between  er  and  F  is  monotonous,  and  log  K2  ~  log  K\ 
will  still  appear  as  a  line  with  slope  of  —1,  and  only  certain 
(usually  slight)  difference  in  the  intercept  compared  with  the 
present  simplified  results. 


3.2.  Effects  of  the  parameters  onflow  distribution 


According  to  Eq.  (41)  and  Fig.  10,  it  is  obvious  that  in  order  to 
reduce  F  and  obtain  more  uniform  flow  distribution,  the  product 
of  K\  and  K2  should  be  reduced.  For  the  cases  with  uniform 
channels  (case  1-0,  2-0  and  3-0),  a  general  expression  for  K1K2 
can  be  given  as  follows. 


KiK2  =  N 


ACDC2 
(Re  f)cLc 


(Re  J%Lh 
AhDh2 


(42) 


It  is  readily  found  that  to  obtain  smaller  F  and  thus  K\  K2 ,  one 
needs  smaller  A,  longer  and  finer  channels  with  larger  aspect 
ratio,  and  shorter  and  coarser  headers  with  smaller  aspect  ratio. 
The  subscript  “h”  here  represents  the  layer  header  when  con¬ 
sidering  the  flow  distribution  among  channels  in  a  layer,  or  the 
stack  header  when  considering  the  flow  among  layers  in  a  stack. 
Since  a  stack  header  is  often  larger  than  a  layer  header,  the  flow 
distribution  among  the  layers  is  usually  more  uniform  than  that 
among  the  channels.  The  variant  um  and  the  fluid  features  (p, 
fi)  have  not  appeared  in  the  above  expression,  and  thus  these 
parameters  cannot  affect  the  flow  distribution  as  reported  by 
other  authors  [5]. 

For  the  channels  consisting  of  two  different  parts  (both 
“cl  +  c2  +  cl”  and  “cl  +  c2”  types,  i.e.,  casel-1,  1-2,  2-1,  2-2, 
3-1  and  3-2),  the  product  of  K\  and  K2  is  expressed  in  Eq.  (43), 
and  the  two  parts  now  both  contribute  to  the  flow  distribution. 
And  according  to  the  expression,  each  part  needs  to  be  fine,  long 
and  of  large  aspect  ratio  to  achieve  a  uniform  flow. 


K\K2  =  N - ^ - y 

(Ref)clLcl/AclD2cl+(Ref)c2Lc2/Ac2D22 

(ReJXLh 

AhD2 


(43) 


Table  1 

The  flow  distribution  in  different  cases  (dimensions  in  mm) 


Case  1-0 

Mi- layer 

ffh-layer 

Mi- layer 

Eq 

Wc 

Eq 

Nc 

KiK2 

Fc 

40 

8 

1 

150 

2 

1 

20 

1.9661 

0.7097 

80 

4 

2 

150 

2 

1 

20 

0.6667 

0.4061 

80 

3 

2.5 

150 

2 

1 

20 

0.6243 

0.3888 

80 

4 

2 

150 

2 

0.5 

20 

0.1033 

0.0886 

80 

4 

2 

150 

1 

1 

20 

0.2080 

0.1641 

Casel-1  or  case  1-2 

Mi- layer 

Mi-layer 

Mi-layer 

Eq\ 

Eq2 

Wc 

Eq\ 

Eq2  Nq 

KxK2 

Fq 

80 

4 

2 

20 

130 

2 

0.5 

1  20 

0.3860 

0.2750 

80 

4 

2 

50 

100 

2 

0.5 

1  20 

0.2366 

0.1852 

80 

4 

2 

75 

75 

2 

0.5 

1  20 

0.1789 

0.1455 

Case2-0 

Mi-  stack 

Wh-  stack 

Eh- stack 

Eq 

Wc 

Eq 

Mayer 

KiK2 

Mayer 

100 

2 

18 

150 

2 

1 

20 

0.1353 

0.1119 

100 

2 

18 

150 

2 

0.5 

20 

0.0210 

0.0187 

100 

2 

78 

150 

2 

1 

20 

0.0297 

0.0264 

Case2-1  or  case2-2 

Mi-  stack 

Mi-stack 

Mi-stack 

Eq  1 

Eq2 

Wc 

Eqi  Ec2  N\  ayei 

KiK2 

F layer 

100 

2 

18 

20 

130 

2 

0.5  1 

20 

0.0783 

0.0671 

100 

2 

18 

50 

100 

2 

0.5  1 

20 

0.0480 

0.0420 

100 

2 

18 

75 

75 

2 

0.5  1 

20 

0.0363 

0.0321 

Case3-0 

Eh-  stack 

Mi-stack 

Eh-  stack 

Eq 

Eq 

Mayer 

KiK2 

Mayer 

100 

8 

8 

150 

2 

1 

20 

0.0104 

0.0094 

100 

4 

4 

150 

2 

1 

20 

0.1669 

0.1351 

100 

4 

4 

150 

2 

0.5 

20 

0.0259 

0.0230 

Case3-1  or  case3-2 

Eh-  stack 

Mi-stack 

Mi-stack 

Eqi 

Eq2 

Wc 

EqI  Eq2  Mayei 

K\K2 

Mayer 

100 

4 

4 

20 

130 

2 

0.5  1 

20 

0.0967 

0.0817 

100 

4 

4 

50 

100 

2 

0.5  1 

20 

0.0593 

0.0515 

100 

4 

4 

75 

75 

2 

0.5  1 

20 

0.0448 

0.0393 
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term  while  retained  the  inertial  term.  However,  with  the  assump¬ 
tion  of  laminar  flow,  Maharudrayya  et  al.  [5]  have  neglected  the 
inertial  term,  retaining  the  frictional  term.  The  resultant  analyti¬ 
cal  results  can  present  some  discrepancy  in  comparison  with  the 
CFD  simulations  without  such  ignorance.  Since  this  term  fur¬ 
ther  reduces  the  pressure  along  the  header,  it  should  increase  the 
flow  nonuniformity.  On  the  other  hand,  this  flow  inertial  move¬ 
ment  might  cause  inverse  flow  at  the  header  end,  causing  certain 
increase  in  mass  flow  rate  in  the  last  several  channels  or  layers, 
even  breaking  the  monotonic  distribution.  Another  assumption 
is  the  ignorance  of  local  pressure  losses,  and  even  the  layer- 
header  losses  within  the  stack  formulations.  Such  losses  might 
be  favorable  for  the  corresponding  flow  uniformity  among  the 
channels  or  layers  by  reducing  K\.  These  competitive  factors 
contribute  to  the  difference  situations  when  comparing  the  CFD 
(without  the  above-mentioned  assumptions)  and  the  analytical 
results,  according  to  the  exact  configurations  and  even  operating 
parameters. 

Maharudrayya  et  al.  [5]  have  provided  certain  CFD  simu¬ 
lations  to  verify  the  validation  of  their  analytical  formulations. 
Although  their  results  exhibited  good  agreements,  the  differ¬ 
ences  can  still  be  observed  to  some  degree  (Figs.  6  and  7  in  Ref. 
[5]).  In  the  first  several  channels,  the  CFD  derived  flow  rates 
are  somewhat  higher  that  the  analytical  ones  mainly  due  to  the 
enhanced  pressure  loss  by  the  inertial  term,  while  in  the  last  sev¬ 
eral  channels,  the  CFD  flow  rates  also  give  higher  values  and 
somewhat  change  the  monotonic  decreasing  of  the  flow  distri¬ 
bution  along  the  feed  header.  Based  on  certain  assumptions,  the 
analytical  formulation  shows  that  in  the  U-type  configuration, 
the  flow  distribution  is  independent  of  the  Reynolds  number. 
However,  slight  difference  appear  in  the  distribution  plots  with 
different  Reynolds  number,  and  higher  Reynolds  number  gives 
more  obvious  effect  from  the  inertial  term  (Fig.  7(b)  in  Ref. 
[5]). 

Figs.  11  and  12  illustrate  the  comparison  between  the 
analytical  and  CFD  (FLUENT™  6.0)  results  in  the  present 
case  1-0,  easel -1  and  case  1-2,  with  common  header  length 


Channels 

Fig.  11.  Comparison  of  flow  distribution  obtained  from  analytical  and  CFD 
solutions  in  case  1-0. 


Channels 

Fig.  12.  Comparison  of  flow  distribution  obtained  from  analytical  and  CFD 
solutions  in  easel -1  and  easel -2. 


(Lh-iayer  =  80  mm)  and  channel  number  (Nc  =  20).  Other  val¬ 
ues  are  Wh-iayer  =  4  mm,  //h-iayer  =  2  mm,  Wc  =  WY. iayer  =  2  mm. 
Additionally,  Lc  =  150  mm  and  Hc  =  1mm  in  case  1-0,  while 
in  casel-1  and  casel-2,  Lci=50mm  (partitioned  into  two 
25  mm  in  casel-1),  Hc\  =  0.5  mm,  Lc 2  =  100  mm,  Hc\  =  1  mm. 
The  Reynolds  number  was  fixed  at  137  for  all  the  cases.  In 
the  present  three-dimensional  CFD  simulations,  a  laminar  flow 
was  assumed  with  constant  p  and  /x,  and  mass  conservation 
(no  gas  consumption).  The  conservation  equations  for  mass  and 
momentum  (the  Navier-Stokes  equations)  were  solved  with  a 
segregated  solver,  and  the  local  losses  and  the  inertial  term  were 
not  ignored.  The  analytical  mass  flow  rates  were  taken  from 
the  channel  centers,  while  those  in  the  CFD  simulations  were 
obtained  via  integrating  over  the  channel  cross-sections.  The 
overall  trends  in  the  two  results  are  in  agreement,  indicating  the 
certain  validation  of  the  analytical  formulations.  The  Fc  values  in 
case  1-0  are  0.4061  and  0.5045  for  the  analytical  and  CFD  results, 
respectively,  and  the  last  channel  (no.  20)  in  CFD  does  not  offer 
the  minimum  mass  flow  rate.  This  discrepancy  is  consistent 
with  the  reports  in  Ref.  [5]  with  different  layer  sizes.  Fig.  12 
gives  the  results  for  casel-1  and  casel-2.  The  analytical  Fc  is 
0.1852,  while  CFD  simulations  give  Fc  =  0.3160  and  0.3323 
for  casel-1  and  casel-2,  respectively.  Although  the  deviation 
is  obvious,  the  size  change  in  channels  is  found  to  be  effective 
in  improving  the  flow  uniformity.  And  the  casel-1  results  from 
CFD  show  better  uniformity  than  casel-2,  indicating  the  pres¬ 
sure  loss  difference  in  the  two  configurations  (the  total  pressure 
drops  from  the  CFD  simulations  are  20.670  and  20.516  Pa  for 
casel-1  and  casel-2,  respectively,  and  for  comparison,  the  CFD 
value  in  casel-0  is  12.024 Pa),  which  has  been  ignored  in  the 
analytical  formulations.  With  the  reduction  in  header  length  and 
channel  number,  the  deviation  due  to  the  inertial  term  can  be 
expected  to  be  more  apparent.  And  for  stack  cases,  the  reduc¬ 
tion  in  header  length  and  layer  number  might  lead  to  more  severe 
discrepancy  between  the  results  with  and  without  considering 
the  inertial  term,  due  to  the  large  cross-sections  of  the  stack 
header. 
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4.  Conclusion 

The  present  work  extends  the  analytical  formulations  for 
uniform  channels  in  a  U-type  layer  to  more  complex  channel 
structures  and  U-type  stacks  based  on  certain  assumptions.  The 
analytical  results  indicate  simply  that  the  flow  distribution  can 
be  improved  by  reducing  the  pressure  loss  in  the  header  (via 
shortening  the  header,  expanding  the  hydraulic  diameter  of  the 
header,  reducing  the  channel  or  layer  number,  etc.),  or  increas¬ 
ing  the  pressure  loss  in  the  channels  or  layers  (via  lengthening 
the  channels,  reducing  the  hydraulic  diameter  of  the  channels 
and  so  on).  According  to  the  formulations,  the  flow-distribution 
uniformity  among  the  layers  in  a  U-type  stack  is  easier  to  obtain 
compared  with  that  among  the  channels  in  a  U-type  layer  due 
to  the  larger  size  of  stack  headers  usually  adopted.  Modifying 
the  channel  shape  with  narrow  parts  can  effectively  improve  the 
flow  uniformity,  according  to  both  the  analytical  formulations 
and  CFD  simulations.  The  discrepancy  between  the  analytical 
and  CFD  results  might  be  ascribed  to  whether  the  inertial  term 
and  local  pressure  losses  have  been  ignored  or  not. 
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